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Abstract 

We reconsider thin-disc global asymptotics for kinematic, axisymmetric 
mean-field dynamos with vacuum boundary conditions. Non-local terms aris- 
ing from a small but finite radial field component at the disc surface are con- 
sistently taken into account for quadrupole modes. As in earlier approaches, 
the solution splits into a local part describing the field distribution along the 
vertical direction and a radial part describing the radial (global) variation 
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of the eigenfunction. However, the radial part of the eigenfunction is now 
governed by an integro-differential equation whose kernel has a weak (log- 
arithmic) singularity. The integral term arises from non-local interactions 
of magnetic fields at different radii through vacuum outside the disc. The 
non-local interaction can have a stronger effect on the solution than the local 
radial diffusion in a thin disc, however the effect of the integral term is still 
qualitatively similar to magnetic diffusion. 

KEY WORDS: Mean-field dynamos, Thin-disc asymptotics, Boundary con- 
ditions, Galactic magnetic fields 



1 Introduction 

In the standard thin-disc asymptotic approach to mean-field dynamos, the dynamo 
equation 

f)F} 

— = Vx (V X B + aB-(3Vx B) , (1) 
at 

splits into a local part containing only derivatives with respect to the vertical coordi- 
nate z and the horizontal part containing derivatives with respect to the cylindrical 
radius r and the azimuthal angle 0. Here B is the mean magnetic field, V is the 
mean velocity, a is the coefficient describing the alpha-effect, and (3 is the turbulent 
magnetic diffusivity. We use cylindrical polar coordinates (r, 0, z) and consider ax- 
isymmetric solutions of the kinematic dynamo equation, d/d(f) = 0. In terms of the 
azimuthal components of the magnetic field B and the vector potential A, Eq. (1) 
reduces to 

— DG— — A 2 — (-— vb) (2) 

dt dz dz 2 dr \r dr j 

d _A aB ?A a 2 — (-— ta) (3) 

dt dz 2 dr \r dr J ' 

where D = aoGh 3 /(3 2 is the dynamo number, G = rdVt/dr (with Q the angular ve- 
locity), and we have used dimensionless variables defined by the disc semi-thickness 
h for z, characteristic disc radius R for r, magnetic diffusion time h 2 / f3 for t, and 
the ratio of units for B and A is aoh 2 /j3, with a a representative value of the 
a-coefficient (see, e.g., Ruzmaikin et al, 1988 for details). We have neglected the 
term with dQ/dz in Eq. (2) assuming that \dQ/dz\ \dQ/dr\ under typical con- 
ditions in astrophysical discs. In terms of the dimensionless variables, we consider 
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the region \z\ < 1, 0<r<oo. Astrophysical discs have no well defined radial 
boundary, and the dynamo active region has a finite size because the intensity of 
the dynamo action decreases at large r; therefore, the characteristic radius of the 
dynamo region R, introduced below, is finite. The disc boundary is not sharp at its 
surface as well, but the turbulent magnetic diffusivity rapidly increases with height 
(at least, in spiral galaxies), so that vacuum boundary conditions posed at a finite 
height are a good approximation (Moss et al, 1998). 

Equations (2) and (3) are written for the aw-dynamo, but generalization to a 2 u- 
dynamos is straightforward as long as the radial scale of the solution in not affected 
much. Our basic result, Eq. (32), remains valid for both aoo- and a 2 a;-dynamos, but 
7o(r) and r](r) defined below will follow from different local equations. 

An asymptotic solution of Eqs. (2) and (3) is represented by Ruzmaikin et al. 
(1988) in the form 




exp (Ft) 



where T is the growth rate of the field (r is known to be real for the dominant modes 
in a thin disc if D < 0), (b, a) is the local solution (suitably normalized) depending 
on r only parametrically, and Q(r) is the amplitude of the global solution. The 
asymptotic parameter is the disc aspect ratio 

A = 4 - lO-MCT 1 . 
R 

Equation (1) should be supplemented by appropriate boundary conditions (dis- 
cussed in detail below), and then the lowest-order approximation in A yields a 
boundary value problem for (b,a), whereas the solvability condition for the first- 
order equations leads to an equation for Q and T. The procedure for the derivation 
of the asymptotic equations is reconsidered in detail below. 

The boundary conditions traditionally used in conjunction with Eqs. (2) and 
(3) are the so-called vacuum boundary conditions. If there are no electric currents 
outside the disc, then VxB = and, for axisymmetric solutions, (see, e.g., Zeldovich 
et al, 1983, p. 151) 

B = for \z\ > 1 . 
Further, V 2 A = outside the disc, so that 



dz 2 dr 
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In order to impose boundary conditions at the disc surface, we recall that the eigen- 
functions of Eq. (4) are of the form A oc exp (—k\\z\)Ji(kr) with some radial wave 
number k. Since we require that the scale of the radial variation of A is of order 
unity, k = 0(1), we obtain 



It was noted by Soward (1992a,b) that the term of order A in Eq. (5) must be 
retained if higher-order terms of the asymptotic solution of Eqs. (2) and (3) are 
considered. Our goal here is to include these terms in a consistent manner into a 
regular asymptotic scheme. The physical meaning of these corrections is as follows. 
Since —dA/dz = B r , the boundary condition (6), together with B = 0, restricts 
the magnetic field to be purely vertical outside the disc, similarly to an external 
field of an infinite homogeneous slab. Then magnetic fields at different radii can be 
connected only via magnetic lines passing through the disc where B r ^ 0. However, 
the boundary condition (5) implies that the field on and above the disc surface has 
a weak but finite radial component, and magnetic lines can thus connect different 
regions in the disc through the vacuum. This leads to non-local effects arising 
in higher orders of asymptotic expansion, and a description of these effects is our 
subject here. 

A treatment of the non-local asymptotics for mean-field dynamos in a thin disc 
can be found in Soward (1978, 1992a,b, 2000) where both steady and oscillatory 
dipolar and quadrupolar modes are considered, but only for marginally stable solu- 
tions and in a local, Cartesian geometry. Here we consider exponentially growing 
(or decaying), non-oscillatory, quadrupolar modes in cylindrical geometry, ready for 
applications to astrophysical discs, especially those in spiral galaxies. Another new 
element in our approach is that we do not expand the coefficients of the dynamo 
equations into power series in r, and so we allow for arbitrary radial variation in 
the dynamo coefficients. In this paper we present the derivation of the asymptotic 
equations. Their solutions will be discussed elsewhere (Willis et al, 2000). 



8A 

dz 



kXA = O(A) for z>l. 



(5) 



To the lowest order in A, the resulting boundary condition is 




(6) 
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2 Vacuum boundary conditions for a thin disc 



In this section we derive boundary conditions for a thin disc surrounded by vacuum 
in a form suitable for asymptotic analysis. For this purpose we first perform the 
Hankel transform of Eq. (4) to obtain: 

g-A^4 = 0, (7) 

where ^ 

A(z,k)= A(z,r)J 1 ( y kr)rdr (8) 
Jo 

is the Hankel transform of A(z,r) with respect to r. Equation (7) has a simple 
solution decaying at infinity, A(z,k) = A (k) exp (—k\\z\). However, this solution 
contains an arbitrary function of k, A (k), so it cannot be used to formulate a closed 
boundary condition at \z\ = 1. Instead, we note that, for z > 0, solution of Eq. (7) 
decaying at z — > oo satisfies 

8 A 

— + XkA = for z > 1 , 
az 

which is useful to compare with Eq. (5). We apply the inverse Hankel transform to 
this equation and set z — 1 to obtain 



dA 

dz 



roc ^ 

+ A / A(l,k)J 1 (kr)k 2 dk = , (9) 
Jo 



and the plus sign must be replaced by minus for z — — 1. 

We could now substitute Eq. (8) into Eq. (9) in order to obtain a closed boundary 
condition for A at z = 1 in an integro-differential form, but the resulting operator 
contains a strongly divergent kernel / °° Ji(kr)Ji(kr')k 2 dk. Therefore, we develop a 
slightly more elaborate procedure devised to isolate the singularity in the integral 
kernel. Firstly, we improve the convergence of the kernel by eliminating the factor 
k 2 in it. This can be done by rewriting the integral operator in terms of d 2 A/dr 2 . 
For this purpose, we use the general relation (Sneddon, 1951, Sect. 10) 

d (1 d 



^ poo y i 1 o \ 
k 2 A = - — — ( — — rA jAkr)rdr , 
Jo or \r or I 
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written at \z\ = 1, to reduce Eq. (9) to the form 



OA 

dz 



-- r rt> {? t> r ' A ) w{r > r,) dr>= ° at z = 1 ' (io) 



where 



poo 

W(r,r') = rr' / JAkriJAkr') dk ; (11) 
Jo 

note that we have isolated a factor rr' in W(r, r') to ensure that it does not diverge 
at r, r' — > (see below). 

We can prove that W(r,r') is square integrable and has only a weak singularity 
at r — r' = 0. For this purpose we explicitly isolate the singularity of W(r, r') which 
occurs at r — r' = 0. We shall consider consecutively the cases prone to a singular 
behaviour. 

We start with r,r' ^ when a singularity at r = r' can be expected. We 
introduce a new variable a = kr for r' > r or a = kr' for r' < r in Eq. (11) and use 
Eq. (Al) of Appendix A to see that 

W(r, r') = --(rr') 1/2 In |r - r'\ + . . . ; (12) 

7T 

here and henceforth dots stand for a nonsingular part. In a similar way the case 
when r, r' — > can be shown to lead to the singularity of Eq. (12) because terms 
like (rr') 1 / 2 In r are not singular. 

Next consider r = const ^ 0, r' — > 0. The new variable is now cr = At', so that 

W(r,r') = (rr') 1/2 In |r - r'\ + -(rr') 1/2 lnr' + . . . . 

7T 7T 

Both logarithmic terms are finite due to the factor rr' in W(r,r'), so in this case 
the kernel is free of singularities. As r and r' enter W in a symmetrical manner, we 
conclude that there is no singularity at r — > 0, r' = const 7^ 0, as well. 

Altogether, Eq. (12) shows that the integral kernel has an integrable (logarith- 
mic) singularity, and thus it can be treated as a finite kernel (see, e.g., Kolmogorov 
and Fomin, 1957). As a result Eq. (10) provides a viable exact boundary condition 
for a slab surrounded by vacuum. 



6 



3 Representation of the Asymptotic Solution 

In this section we consider asymptotic solutions of Eqs. (2) and (3) which are free of 
boundary layers. Such layers may arise at the disc surface and also at certain radii 
(e.g., contrast structures — Belyanin et al, 1994). 

/ Biz r) \ 

With the definition Y — I .) ' ! , the eigenvalue problem for Eqs. (2) and 

\ A{z,r) J 

(3) can be conveniently rewritten in the form 

FY = C Z Y + X 2 C r Y , (13) 

with the boundary conditions 

V zfi Y = at z = ; V Z>1 Y - ^KY = at z = 1 , (14) 

where the following operators have been introduced: 

_ / dydz* -DGd/dz \ d (Id \ 

Lz ~{ a d 2 /dz 2 )> LrA ~ dr [r dr ) ' 

V _(d/dz 0\ _(1 \ 



and 



x 2 ) y J™ dr'W{r,r')C r ,{x 2 ) 



The operators V z> i and /C are defined at z — 1 and "P Zj0 at z = 0. 

Consider asymptotic solutions of Eqs. (2) and (3) for A <C 1. We denote the 
asymptotic parameter e, a function of A to be determined below, and represent the 
asymptotic solution in the form 

Y(z,r) = Q(e- 1 r)Y (z;r) + e 2 Y 1 (z,e-' 1 r;r) + ... , (16) 

with unknown functions 

r>(,-K\ v - ( 6 ^ ;r ) ^ v - ( B i( z i € ~ lr ; r ) \ 

g(£ r) ' Y °-[a(z;r)h Yl ~ [ A ± (z, e^r; r) J " 
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Our motivation for the particular choice of this form will immediately become clear. 
We only note here that Eq. (16) features two types of dependence on r. One of them 
is connected with the presence of the radial derivatives in Eq. (13), this involves the 
argument e _1 r. Another kind of r-dependence arises from the radial variation of 
the coefficients of the local operator, C z ; this is a slow parametric variation of the 
solution with r unrelated to any differential operator in r. The asymptotic expansion 
of the eigenvalue T will emerge later. 

In order to preserve a similarity to earlier asymptotics developed for galactic 
dynamos, we do not introduce scaled variables, e.g., e _1 r. Therefore, some of our 
asymptotic equations, e.g., Eq. (32), will formally contain terms with e or A, but all 
terms in such equations will always be of the same order of magnitude in A. 

3.1 Zeroth-order equations 

Substituting Eq. (16) into Eqs. (13) and (14) and combining terms of equal orders 
in e we obtain in the leading order in A 

(TY o -C g Y o )Q = 0, (17) 

QP Z , Y = , QV zA Y = , (18) 

where we have deliberately retained the factor Q, a function of e _1 r. If Q were 
cancelled, Eq. (17) would become intrinsically contradictory as T must be indepen- 
dent of r whereas the coefficients of C z depend on r. To overcome this difficulty, 
we consider separately two radial ranges. At those r where Q is small, Eqs. (17) 
and (18) are satisfied, to the required accuracy of 0(e 2 ), just because Q is small, of 
order e 2 or less. At radii where Q = 0(1) we replace Eqs. (17) and (18) with the 
following equations: 

7o(r)Yo-£,Yo = 0, (19) 

V Z , Y = , V zA Y = , (20) 

where 7o(?") is the local growth rate (more precisely, its zeroth-order term), a function 
of radius. With this definition of 7o(r), we have 

r- 7o (r) = 0(e 2 ) if Q^r) = O(l) . (21) 

Outside the radial range where Q(e^ 1 r) = 0(1), 7o(r) still can be defined as the 
leading eigenvalue of Eqs. (19) and (20). 
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Equation (21) in fact restricts the radial range Ar where Q = 0(1) to be Ar = 
0(e) provided 70 (r) can be approximated by a quadratic function of r near its 
maximum, 70 (r + Ar) ~ 70 (r) + CAr 2 with a certain constant C. The radial 
variation of 70 (r) is controlled by that of a and G, the coefficients of C z . The weaker 
is the radial variation of a and G, the wider is the radial range where Eq. (21) can 
be satisfied. 

Earlier experience with solutions of Eqs. (19) and (20) (Ruzmaikin et ai, 1988) 
shows that the leading local eigenvalue is real for D < and, furthermore, this 
eigenvalue is isolated, i.e., there are no other positive eigenvalues for moderate values 
of \D\ (in other words, the difference of the eigenvalues is of order unity). We restrict 
our attention to the mode with the largest real eigenvalue 7o(r). 

Since the boundary value problem (19) and (20) has been extensively studied 
earlier (e.g., Ruzmaikin et ai, 1988 and references therein), we do not discuss it 
further, but assume that its solution, Y (z;r) and 70 (r), is given. 

3.2 First-order equations 

Having isolated zeroth-order terms, we can represent the remaining terms in Eq. (13) 
as 

TY 1 = C Z Y 1 - T ~] oir) Y Q + ^Y C r Q + . . . , (22) 
e z e z 

where dots denote terms of order e 2 and higher. We temporarily retain the term of 
the order A 2 /e 2 as the connection of e and A still has not been determined; we shall 
show that this term must be omitted in the first-order equations. 
Similarly the terms remaining in the boundary conditions are 

V Z , Y 1 + ... = 0, 

V^Y, = ^-IC(QY ) + ... . 

e z r 

We determine e from the requirement that the two terms in the second boundary 
condition are of equal orders of magnitude in A. As shown in Eq,. (16), the eigen- 
function is expected to be of order unity in a region whose radial extent Ar is of 
order e, so that dQ/dr ~ Q/Ar ~ Q/e. This is true for the leading eigenf unctions 
Q that do not rapidly oscillate with r. Since C r A = 0(Ar~ 2 ) = 0(e~ 2 ), we obtain 
IC(QY ) ~ e~ 2 Ar ~ e _1 (see Appendix B). Thus, we require that 




9 



or 

e = A 1/3 . 

Then Eq. (21) implies that 

r-[7o(r)]| max = 0(A 2 / 3 ). 

Now we can see that A/e = o(l), so that A 2 e - 2 £ r (Q) = 0(A 2 / 3 ) and the desired 
first-order boundary value problem based on Eq. (22) reduces to 

r*i = c z y, - r y / ° 3 (r V g , (23) 

V Z , Y 1 = , V Z , 1 Y 1 = —JC(QY ) . (24) 

r 

However, this form of the first-order boundary value problem is inconvenient for 
further analysis. We rewrite it to obtain a boundary value problem with homoge- 
neous boundary conditions. Let us represent the solution of Eqs. (23) and (24) in 
the form 

Y 1 = X 1 + X 2: 
where Xi solves the following boundary value problem: 

c z x, = , 

V Z , X 1 = , V Z , 1 X 1 = — /C(QF ) . 

r 

Then X 2 is the solution of the following boundary value problem 

TX 2 - C Z X 2 = 7o( ^~ r y Q - YX l , (25) 

V z , X 2 = 0, V z , l X 2 = 0. (26) 

Here Xi is in fact a first-order correction to the eigenfunction Y of Eqs. (17) 
and (18) arising from the term with A 1 / 3 in the boundary condition (24) at z — 
1. Note that \X ± \ = 0{Q). The boundary value problem for X x can be solved 
straightforwardly upon solving Eq. (32) for q(r) = Q(r)a'(l,r). 

It is useful to represent Y\ as a sum of the two terms, Xi and X 2 , because this 
results in a universal form of the left-hand sides in the local boundary value problems 
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in all asymptotic orders. An alternative approach, without this representation, is 
discussed by Soward (2000). 

We expect that X 2 has the same functional form as the zeroth-order term in 
the asymptotic expansion (16), so we take X 2 (z;r) = Qi(\~ 1 / 3 r)X 2i0 (z;r). Then 
Eqs. (25) and (26) reduce to a form similar to Eqs. (17) and (18): 



7o(r) — r 

T^2,o — (£zX 2l o)Qi = — 7575 — Y Q — TXi , 



A 2 / 3 

QiP Z fiX 2 fi = , QiP z ,\X 2 fi 







(27) 



f28) 



where, as above, we retain the factor Qi in the boundary conditions and consider 
again two radial ranges. At those radii where V — 7o(r) = 0(A 2 / 3 ) and Q = 0(1), 
we expect that Qi = 0(1). In this radial range we obtain the following first-order 
equations: 

7o(r) - r. 



rx 2 ,o-£ 2 x 2 , = ,uv A 2 /3 : y oQ - rxi Qr 1 , (29) 

V g ,oX 2fi = , V g ,iX 2fi = . (30) 

These equations have the same operators on the left-hand sides as Eqs. (19) and 
(20). As for the radial range where Qi = o(l), the form of Q\ is unimportant there 
and we do not discuss it further. Equations (29) and (30) represent the desired form 
of the first-order boundary value problem which we use below to obtain an equation 
for Q. 



4 The radial dynamo equation 

In this section we derive a closed equation for the radial part of the zeroth-order 
eigenfunction, Q, which follows from the solvability condition for the first-order 
boundary value problem (29) and (30). By Fredholm's solvability condition, the 
r.h.s. of Eq. (29) must be orthogonal to the eigenvector X^ of the corresponding 
adjoint homogeneous problem (this problem is formulated and discussed in Ap- 
pendix C), i.e., 

where (F, G) — JqF-G dz denotes the scalar product in the space of continuous 
vector functions on the interval < z < 1. This yields the following equation for Q 



r -7o(r) 

A2/3 



Y Q + TX 1 
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and T: 

T ~ x y} r) Q(Y , xt) + r(x 1 , xt) = o . (3i) 

We recall that T — 7o(r) = 0(A 2 / 3 ) and \Xi\ = 0(1) in the radial range where 
Q = 0(1) and \Xi\ = o(l) outside this range. Therefore, to the accuracy 0(A 2 / 3 ) 
we can replace T by 70 (r) to obtain 

r(X!,Xt) = (X 1)7 o(r)Xt) 

= (^i,4^> 

\l/3 

= -a f (l,r)--/C(aQ) , 

where, in order to obtain the last equality, we performed two integrations by parts 
and took into account the boundary conditions (20) and (C3). 
Now Eq. (31) can be rewritten as 



TQ = 7o (r)Q + *° t(1; J° t , r dr> W(r, r')L, (Qa(l, r'] 



This is the desired equation for magnetic field amplitude Q and global growth rate 
T. We stress that all terms in this equation are of the same order in A because 
dQ/dr = 0(\- 1 / 3 ). 

To reduce this equation to a final form, we introduce 

q(r)=Q(r)a'(l,r) , 

to obtain 

A f°° 

Tq = 7o (r)g + -r](r) / dr' W{r, r')C r , (q) , (32) 
r jo 

where 



a(l,r)at(l,r) 
" (r) = { Y ,X<) 
is a function of radius which is obtained from the zeroth-order equations, 

_ / b(z;r) \ ( b\z;r) \ 

° [a(z-r) J ' A \ct(z;r) J ' 

It is important to note that both q(r) and i](r) are independent of the normalization 
of the local solutions. 
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5 Discussion 



It is useful to compare the main result of this paper, Eq. (32), with the corresponding 
equation resulting from using the boundary condition (6) in all orders (see, e.g., 
Ruzmaikin et al, 1988): 

rg = 7o (r)Q + \ 2 C r Q . (33) 

The last term on the right-hand side describes magnetic diffusion in the radial di- 
rection. The new feature of the theory discussed here is that this term is replaced 
by an integro-differential operator. Numerical analysis of Eq. (32) presented else- 
where (Willis et al., 2000) shows that the integral term results in enhanced magnetic 
diffusion, so that Eqs. (32) and (33) have similar physical significance. This is under- 
standable since the kernel W(r, r') has a singularity. If this singularity were extreme, 
say, W(r,r f ) ~ S(r — r'), the functional form of Eq. (33) would be recovered pre- 
cisely, albeit with additional factor i](r) which can hardly be of major importance. 
However, the actual kernel has a wider radial profile (and weaker singularity), so 
the magnitude of this term is enhanced by non-local coupling, but its nature still 
remains basically diffusive. 

The integral, non-local character of magnetic diffusion affects the nature of the 
asymptotic solution for Q(r). From Eq. (33), the radial width of the eigenfunction 
is of the order of e = A 1 / 2 and T — 7 = 0(A) (Ruzmaikin et al, 1988), whereas 
Eq. (32) corresponds to e = A 1 / 3 and T — 70 = 0(A 2 / 3 ). We should note, however, 
that the difference is hardly important for most applications to astrophysical discs 
where, for A ~ 10~ 2 , the difference amounts to a factor of A 1 / 6 w 1/2. A wider radial 
distribution of the magnetic field resulting from Eq. (32) (Soward, 2000; Willis et 
al, 2000) can help in accelerating the amplification of magnetic fields in the outer 
regions of young galaxies along the lines discussed by Moss et al. (1998). 

Using heuristic arguments, Poezd et al. (1993) generalized Eq. (33) to nonlinear 
regimes arising from alpha-quenching. The resulting nonlinear equation has T re- 
placed by d/dt and 7o(r) by 7o(r)[l — Q 2 / 'B%(r)\, where B (r) defines the dynamo 
saturation level appearing in the alpha- quenching model. It is plausible that a sim- 
ilar generalization is possible with Eq. (32). Another problem is to generalize the 
theory to nonaxisymmetric solutions where the boundary conditions should take a 
more complicated form. 
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Appendix A. The singularity of J Q J\{k)J\(k(j) dk 

In this section we isolate the singularity of the integral 

1(a) = / J 1 (k)J 1 (ka)dk , 
J o 

where a is a constant. It is a special case of the Weber and Schafheitlin integral 
discussed at length by Watson (1922, Sect. 13.4), which can expressed as 

I{a) = \a 2 F l (\,\,2,a"). 

We can restrict ourselves to the case a > 1 because Eq. (11) reduces to the above 
form with a — r/r' for r > r' and a = r'/r for r' > r. 

Asymptotics of this integral for a — 1 <C 1 can be obtained by expressing 1(a) 
in terms of the complete elliptic integral of the first kind, which yields 1(a) ~ 
-it- 1 In 1 1 — cr | (see Willis et al, 2000 for details). 

This asymptotic form can also be obtained from the following arguments which 
make clear the nature of the singularity at a — > 1. Introduce 



A(x) = Ji(x) - ^— j cos (x - f tt) , 
so that A(x) = 0(x~ 1 ) for x > 1. Then 

rl roo 

1(a) = J J 1 (k)J 1 (ka)dk + J A(k)A(ka)dk 

A(k) cos(ka - |tt) f —J + J A(ka) cos(k - f tt) f — J dk 
+ 2^ cos(/c - |tt) cos(ka - |tt) ^ 1/2 , 
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where the only divergent term is the last one. For a > 1, the latter integral can 
be expressed in terms of the sine integral Si(x) = — J X u" 1 smudu and the cosine 
integral Ci(x) = — u^ 1 cos udu (which has a singular component lnx) as 



/oo 
cos(A; 



dk 1 

fvOcos^a-fTr)-^ = — ^ [-|tt - Si(<7 + 1) - Ci(a - 1) 



where f(a) is a finite function. 
Hence, for a > 1 we have 



/°° J^J^ka) dk = F(a) ^ In (a - 1) , (Al) 



where F(a) is finite. 



Appendix B. The order of magnitude of the inte- 
gral term 

Here we estimate the order of magnitude of 

]C(Qa) = Vw(r,r')£ r ,(Q(e-V)a(z;r')) . 

Since a is a slowly varying function of r', its radial derivatives are of order unity. 
The dominant term in C r iQ arises from d 2 Q/dr' 2 = 0(e~ 2 ). Since the magnitude of 
the integral is controlled by the singularity in W(r,r'), Eq. (12), we have 

!C(Qa) ~ [°° \ n \r - r'&r dr 1 . 
Jo dr' 2 

It may seem that the order of magnitude of this integral is e~ 2 Ar In |Ar| ~ e _1 lne 
because the main contribution to the integral arises from a neigbourhood of r = r' 
whose width is Ar ~ e. 

However, this estimate is wrong because d 2 Q/dr' 2 changes rapidly and even 
changes sign within the e-neighbourhood of r — r'. The correct calculation is as 
follows. We introduce 5 — > and consider 

-s d 2 



In r — r — — - dr ~ / m\r — r \—— dr + / m\r — r \—— dr , 

dr 1 Jo dr' 2 J r +s dr' 2 
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so that now we can integrate by parts to obtain 

~1 r— 5 



K{Qa) 



dQ 
dr' 



In |r — r'l 



+ 



o 



dQ 
dr' 



In \r — r'l 



- v.p. 



r+<5 



dQ 



o dr' r — r' 



dr' . 



Since the In \r — r'\ is symmetric with respect to r, r' and dQ/dr = for r = and 
oo, the logarithmic terms cancel and 

dQ 1 



IC(Qa) = —v.p. 



o <ir' r — r' 



■dr 7 



Appendix C. The adjoint problem 

Consider the homogeneous counterpart of the boundary value problem (29) and 
(30), i.e., 

C z X- lo X = 0, (CI) 

V z , X = 0, V z , 1 X = 0. (C2) 

This boundary value problem coincides with the zeroth-order problem (19) and (20). 
Our aim here is to formulate the adjoint problem and to discuss its properties. 
We introduce the eigenvectors of the original and adjoint problems, respectively: 

\a(z;r)J \a^(z;r)J 
By the definition of the adjoint operator, we have 

(C Z X, X f ) = f 1 U-DGa' + b"p + (ab + a")a*] dz 
Jo 1 J 

= {—DGab + b'b - bb V ) |* + (a'a f - aa V ) 

+ f [b (tf" + aa*) + a (patf + cfj 

Therefore, the problem adjoint to (19), (20) is 

7o(r)fet = b^ ' + aa) , 
7 t( r ) a t = DGtf' + a V \ 

a t(0) = , 6 t; (0) = , 6 f (l) = , a t; (l) = . (C3) 
According to Fredholm's theory, the adjoint problem has the same spectrum as the 
problem (19) and (20), i.e., 7oH r ) = 7o- 



+ 



dz 
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